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Abstract. We consider a gas of point particles moving on the one-dimensional line 
with a hard-core inter-particle interaction that prevents particle crossings — this is 
usually referred to as single-file motion. The individual particle dynamics can be 
arbitrary and they only interact when they meet. Starting from initial conditions 
such that particles are uniformly distributed, we observe the displacement of a tagged 
particle at time t, with respect to the initial position of another tagged particle, such 
that their tags differ by r. For r = 0, this is the usual well studied problem of 
the tagged particle motion. Using a mapping to a non-interacting particle system 
we compute the exact probability distribution function for the two-tagged particle 
displacement, for general single particle dynamics. As by-products, we compute the 
large deviation function, various cumulants and, for the case of Hamiltonian dynamics, 
the two-particle velocity auto-correlation function. 
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1. Introduction 

Starting with the pioneering work of Jepsen pQ and Harris [2], the study of tagged 
particle motion has been an area of very active research. Much of the earlier studies 
focused [SlIilElISlIZlIHlinilinillllllSlIISlIIlllISlIinillZlIIHlIiniEn] on the statistics 
of the typical displacement of the tagged particle, in particular on the mean square 
displacement in various single-hle systems. One of the most interesting result is that, 
for a Hamiltonian one-dimensional gas of hard particles which move ballistically between 
elastic collisions, a tagged particle moves diffusively, thus the mean square displacement 
of the tagged particle, in a time duration t increases as (Xf) ~ t. On the other 
hand, for a gas of Brownian particles with hard-core interactions, the tagged particle 
motion is sub-diffusive, with (Xf) oc The typical fluctuations of the tagged particle 
displacement is described by a Gaussian distribution with the above variance. Recently 
there has been interest in studying the probability of atypical fluctuations of the tagged 
particle displacement [2T[ [221 ESI ES]- Some of the results for these simple classical 
interacting particle models have been useful in obtaining results for one-dimensional 
quantum systems [231 ES] • 

A number of different theoretical approaches have been used to study the probability 
distribution of the tagged particle displacement. These include the original ideas of 
Jepsen and Harris of mapping to a non-interacting system, exact solution of multi¬ 
particle Fokker-Planck equation with reflecting boundary condition between neighboring 
particles, and the recently developed approach of Macroscopic fluctuation theory. In our 
recent works [I61E2], we have shown a simpler way (as opposed to earlier approaches in 
[HElll]) of using the non-interacting picture to computing tagged particle statistics. In 
the present work, we extend this method to study a particular two particle distribution, 
dehned below. 

We consider a hard-point particle system on the inhnite line. The particles are 
distributed uniformly with a hnite density p. For the case of Hamiltonian dynamics (of 
equal mass particles), the initial velocities are taken to be independent and identically 
distributed random variables. The particles move ballistically in between elastic 
binary collisions. During collisions, the two colliding particles merely interchange their 
velocities. As a result, any set of particle trajectories of the interacting system can 
be constructed from the set of trajectories of a non-interacting system (where the 
trajectories pass through each other) — by exchanging the identities (tags) of the 
particles at crossing. For the case of Brownian particles with hard-point interactions, 
Harris [2] defined the interacting particle problem by starting with the non-interacting 
trajectories and exchanging particle identities whenever two trajectories cross. This 
dehnition is equivalent to enforcing reflecting boundary conditions between nearest 
neighbor pairs in the full multi-particle propagator (thus each particle acts like a hard 
reflecting wall for its nearest neighbors). 

One can generalize this collision rule to other cases where the individual particle 
dynamics is neither Hamiltonian or Brownian, for example Levy walks or fractional 
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(Non-interacting) (Interacting) 

Figure 1. (Color online) An interacting hard-point particles system can be 
constructed from a non-interacting system by exchanging tags (colors) when two 
trajectories cross. 


Brownian motion. Hence in general we define the interacting problem as follows: 
start with the non-interacting trajectories and interchange particle labels whenever two 
trajectories cross. 

For this general interacting particle system, let us tag two particles with tag indices 
0 and r, such that there are r — 1 particles in between them. Let the position of 
these particles be xfit) and xfit) respectively. Here we consider the displacement 
Xt = xfit) — a;o(0) and compute it’s statistics. 

2. Main steps of the calculation 

Initially, we consider 2N + 1 particles, independently and uniformly distributed in the 
interval [—L, L] and evolve them on the inhnite one dimensional line. Since during a 
collision each particle acts as a reflecting hard wall for the other and the particles are 
identical, one can effectively treat the system of the interacting hard-point particles as 
non-interacting by exchanging the identities of the particles emerging from collisions 
[see figure and discussion in previous section]. In the non-interacting picture, each 
particle executes an independent motion and the particles pass through each other when 
they ‘collide’. The position of each particle at time t is given independently by a single¬ 
particle propagator of the general form 

G{y, t\x, 0) = — / (-—, (1) 

o-t \ at J 

where f{—w) = f{w) > 0 and {\y — x\)/at = \w\f{w) dw = A is finite. Evidently, 

/(tc) dw = 1. The dependence on time only appears through the characteristic 
displacement at in time t. While for stochastic processes the propagator arises naturally, 
for Hamiltonian systems (where the dynamics is deterministic) it comes from the 
distribution, taken to be of the form v~^f(v/v), from which the initial velocities of 
the particles are chosen independently. In many problems of interest, the propagator 
happens to be Gaussian, i.e., f(x) = , and a/ is the variance. For example, 

for Brownian particles, at = \/2Dt, where D is the diffusion coefficient, while for 
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Figure 2. We mark the position by a; = a;j(0) of the particle at t = 0 and the 
position y = xfit) of the particle at time t. In this example, j = 6 and k = 9. The 
variable Xt = y — x denotes the displacement of the particle at time t with respect 
to the initial position of the particle. 


Hamiltonian dynamics with Gaussian velocity distribution we have at = vt. Similarly 
for fractional Brownian motion, at oc t^, where H is the Hurst exponent. However, our 
analysis is valid for any general propagator. Note that the dependence on time only 
appears through the characteristic displacement at in time t. 

We mark the position x = a:j(0) of the particle at t = 0 and the position 
y = Xk(t) of the particle at time t [see figure [^. Let Xt = y — x, he the difference 
between these two positions. Our goal is to study the statistical properties of this 
random variable Xt, in the thermodynamic limit N ^ oo, L ^ oo while keeping 
N/L = p hxed. In this limit, in the bulk, the statistics of Xt should depend only on the 
difference of the tags r = k — j, rather than the individual tags j and k. Therefore, we 
set j = + 1 (the middle particle) and fc = iV + 1 + r particle counted from the 

middle particle) before taking the thermodynamic limit. In the thermodynamic limit, 
Xt denotes the displacement of a particle at time t, with respect to the initial position 
of another particle such that their tags differ by r. For r = 0, this represents the usual 
problem of tagged particle displacement, studied in [22] . 

In the following we proceed with the calculation using the non-interacting picture 
discussed above. 

2.1. The joint PDF of two particles 

The joint probability density function (PDF) P{x, j,0;y, k,t) of the particle being 
at X at time t = 0, and the particle being at y at time t, can be expressed in terms 
of properties of the non-interacting particles. In the non-interacting picture, there are 
two possibilities: (i) the particle at time t = 0 becomes the kfi^ particle at time t, 
(ii) a second particle becomes the particle at time t [see hgurej^. We need to sum 
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Figure 3. In the non-interacting picture, there are two possibilities: (i) the j**' 
particle at time t = 0 becomes the particle at time t, (ii) a second particle becomes 
the particle at time t. For the case (ii) there are four situations (a) x < x and 
y < y, (b) X > X and y > y, (c) x < x and y > y, and (d) x > x and y < y respectively. 


over these two processes to get, 

= P(i)(a;,j,0;t/,/c,t)+ P(2)(a;,j,0;t/,fc,t) , (2) 

where P(i) and P( 2 ) are the joint PDFs corresponding to the processes (i) and (ii) 
respectively. 

To compnte the contribntion from process (i) we pick one of the non-interacting 
particles at random at time f = 0, mnltiply by the propagator [given in ([^] that it goes 
from (x, 0) to (t/, t), and then mnltiply by the probability that it is the particle at 
t = 0 and the Ifi^ particle at time t. Thus we obtain the corresponding joint PDF as 

(2N + Ij 

P(i){x,j,0]y,k,t) = — ^j^G{y,t\x,0)FiN{x,j,y,k,t), (3) 

where FiN{x,j, y, k, t) is the probability that there are (j — 1) particles to the left of x 
at f = 0 and {k — 1) particles to the left of y at t. 

To compute the contribution from process (ii), we first pick two particles at random 
at time t = 0, and multiply by the propagators that they go from (x, 0) to {y, t) and 
(x, 0) to {y,t) respectively. We then multiply by the probability there are an (j — 1) 

particles on the left of x at time t = 0 and {k — 1) particles to the left of y at t. Finally, 

integrating with respect to x, y, we get the joint PDF corresponding to this process as 

(2iV-b 1)(2A^) ^ ^ ^ 

Pi2){x,j,0;y,k,t) = -- j dx j dy G{y,t\x,0) G{y,t\x,0) F2n{x, j,y, k,x, 

(4) 

where F 2 N{x,j, y, k, x, y, t) is the probability that there are (j — 1) particles on left of x 
at f = 0 and {k — 1) particles on the left of y at time t, given that there is a particle at 
X at time t = 0, and a particle at y at time t. 

To proceed further, we need the expressions for Fi^ and F 2 n- Let p-+{x,y,t) be 
the probability that a particle is to the left of x at f = 0 and to the right of y at time t. 
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Figure 4. A particle goes from the initial position x to the position y in time t. ni 
denotes the number of particles going from the left of x to the left of y, denotes 
the number of particles going from left of x to the right of y, denotes the number 
of particles going from the right of x to the left of j/, and 77,4 denotes the number of 
particles going from the right of x to the right of y, and ni +772 + + = 2-/V. There 

are j — 1 particles on the left of a; at t = 0 and k—1 particles on the left of y at time t. 


Similarly, we define the other three complementary probabilities. Clearly, 


P-+{x,y,t) = (2L) ^ j 

nx 

i dx' j 

poo 

1 dy'G{y'fi\x',0), 

y 

(5a) 

p+-{x,y,t) = (2^)"^ j 

r dx' j 

px 

^ dy'G{y'fi\x',0), 

— 00 

ry 

(5b) 

p—{x,yfi) = (2^)"^ j 

I dx' j 

pL 

1 dy'G{y',t\x',0), 

— 00 

poo 

(5c) 

P++{x,y,t) = (2^)"^ j 

1 dx' j 

dy'G{y' fi\x' 
y 

(5d) 


and P++ + + p _^ + p _= 1. Armed with the above fonr probabilities, we now 

evalnate Fi^r and F 2 N below. 

2.2. Evaluation of Fi]y{x, j,y, k,t) 

In this case, ont of 2N + 1 particles, the particle at the initial time goes from the 
position X to the position y in time t and becomes the particle at the hnal time. The 
remaining 2N particles are independent of each other and the selected particle. Let ni 
be the nnmber of particles going from the left of x to the left of y, n 2 be the nnmber of 
particles going from left of x to the right of y, be the number of particles going from 
the right of x to the left of y, and 774 be the number of particles going from the right of 
X to the right of y [see hgure |^. Clearly ni + 77-2 + ns + 77,4 = 2N. Moreover, since there 
are j — 1 particles on the left of x, clearly, ni + n 2 = j — 1 and ns + n 4 = 2N — {j — 1). 
Similarly, since there are k — 1 particles on the left of y, we have ni + ns = k — 1 
and n 2 + n 4 = 2N — [k — 1). These equalities imply n 4 — ni = 2N + 2 — k — j and 
n3-n2 = k- j. 
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The number of ways of choosing the set {ni, n 2 , ns, n 4 } is given by the multinomial 
coefficient 

(2iV)! 

ni!n2!n3!n4! ’ 

and each possibility occurs with probability p'fi_p^^p']f_p^\. Hence, summing over all 
possible values of {ni, n 2 , ns, n 4 } we get 


-^iv — 


E 


(2iV)! 


ni-\-n2-\-ns-\-n4=2N 


ni!n2!ns!n4! 


plfipffip^-P'lfi S[n4—ni—2N—2+k-\-j] (5[ns—n 2 —/c+j], 

( 6 ) 


where 6[n] is the Kronecker delta function: 5[n] = 1 if n = 0 and 6[n] = 0 for n 7 ^ 0. 
Now, after using the integral representation of the Kronecker delta, 

1 /‘ 6 o+ 27 r 


6[m — n] = — 
27r 


ri)e (where the initial phase ^0 is arbitrary), (7) 


'do 


in the above equation, it immediately follows that 




r‘(/)o+27r 


' <t>0 


d(f) 




^Oo 


ZTT ^ ^ 


where 


H{x,y,9,(f),t) 


p++{x,y,t)e^'f’ + p—{x,y,t)e + p+_{x,y,t)e^^ + p_+{x,y,t)e (9a) 

1 - (1 - COS0) (p++ +p _) +isin 0 (p++ -p _) 

— (1 — cos^) (p+_ + p_+) + i sin 6 ^ (p+_ — p_+). (9b) 


Now, using the fact that 2N is even and the above integral remains unchanged if both 
0 and 6 are shifted by vr simultaneously, the range of the 0 integral has been broken 
into two parts and each of these contributes equally. After appropriately choosing the 
initial phases, this gives. 


FiN{x,y,t) 



[H{x,y,9,(l),t) 


^-i<t>{2N+2-k-j)^-ie(k-j) 


( 10 ) 


2.3. Evaluation of F 2 n{x, j,y, k,x,y,t) 

In this case a particle whose initial tag is different from j, goes from x to y and becomes 
the k^^ particle at time t, while the initial particle goes from a; to y in time t whose 
hnal tag is different from k. To compute F2N, one has to keep track of both these 
particles. Apart from these two particles, let there be ni particles going from the left 
of X to the left of y, n2 particles going from the left of x to the right of y, ns particles 
going from the right of x to the left of y, and n 4 particles going from the right of x 
to the right of y. Since two of the particles are considered separately, the rest can be 
chosen in (2N — l)!/(ni!n 2 !ns!n 4 !) different ways and ni + n 2 + ns + n 4 = 2N — 1. The 
other two constraints among {njj’s are given by n 4 — ni = 2N -\-2 — k — j-\-Xi and 
ns — n 2 = k — j + X2- Unlike the previous case where Xi = X 2 = 0, here their values 
depend on the order of the positions (x, a;) and {y,y). There arises four situations: 
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Figure 5. A particle goes from the initial position x to the position y in time t, while 
at the same time, another particle from its initial position x goes to y. ni denotes the 
number of particles going from the left of x to the left of y, n 2 denotes the number 
of particles going from left of x to the right of y, denotes the number of particles 
going from the right of x to the left of y, and 714 denotes the number of particles going 
from the right of x to the right of y, and 711 + 712 + 773+714 = 27V — 1. There are j — 1 
particles on the left oi x dX t — Q and k — 1 particles on the left of y at time t. 


(a) X < X and y < y, for which Xi = 1 and X 2 = 0, 

(b) X > X and y > y, for which Xi = — 1 and X 2 = 0, 

(c) X < X and y > y, for which Xi = 0 and X 2 = 1, 

(d) X > X and y < y, for which Xi = 0 and X 2 = —1- 

Now following the procedure used to evaluate Fiat, it is easily found that case, one 
has to keep track of the 


F2N{x,y,x,y,t) = I ~ I ^ [H{x,y, 9 ,((),t)Y^ ^e ^(2^+2 ^ dg dg 70xig ^Sx2_ 


de 


' —TTj2 


( 11 ) 


2 . 4 . Exact PDE of the two-tag displacement in the thermodynamic limit 

So far our calculations are exact, valid for any N and L. We now assume both N and L 
to be large and keep only the dominant terms. Finally, we will take the thermodynamic 
limit N —)■ cx), L —>■ cx) while keeping N/L = p hxed. We now set k = j + r and 
i = A^+1. We also change our notation to F(i)(x, iV+1, 0; t/, A^+l+r, t) —)■ F(i)(+, t/, r, t). 
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P{x,y,r,t) = P(i){x,y,r,t) + P^ 2 ){x,y,r,t). (12) 

From ([^ and (10) we get 

\i){x,y,r,t) =-^f{z) [ ^ (13) 


Pr 


(yt 


' —7r/2 


where z = {y — x)/at. Similarly from (|^ and ( |IT| and performing the integration 
X and y, we get 


over 


P^ 2 ){x,y,r,t) =p^ 


dfi r do 


'-7r/2 




X \2Ai{z)A2{z)cos(t)-\-A\{z)e + A\{z)e^^^^^ (14) 

where the fnnctions y4i_2(^) are given by 

poo poo 

^i(^) = / G{x,t\0,0) dx = / f{w)dw, and A 2 {z) = 1 — Ai{z). (15) 


' CTtZ 


Now we explicitly compnte the expressions for p±± nsing ([^. Keeping only the 
dominant terms np to 0(1/L), which snrvive in the limit N —)• oo, L —)• oo while 
keeping N/L = p hxed, we get 


P-+ = 


P+- = 


fi± 

2L 

(yt 


-\ + Q{z) 


+ 


2L 

1 

2 2L 

1 CXt 

P++= n + WT 


^ + Qiz) 


+ 


p— = 


z 


+ 


2 ’ 2L L 2 


+ ■ ■ ■ , 


where x = {y — x)/at, z = {y x)/at, and 

pz poo 

Q{z) = z / f{w)dw-\- / wf{w)dw. 


(16a) 

(16b) 

(16c) 

(16d) 

(17) 


Now, snbstitnting p±± in the expression (9b), for large N, keeping only the most 
dominant terms, one hnds 


lj‘^N _ —2N{l—cos(p) —ipcrtzsm(f) —2pcrtQ{z){l—cos6) ipcrtzsin6 

J. J. o o o o « 

The PDF oi Xt = y — X is given by 


poo noo 


Pa.g(X,,r,t) = j / S(x,-[ii - x\) P{x,y,r,t)dxdy. 


(18) 


(19) 


' —OO J —oo 
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Using (18), in (12)-(14), and making a change of variables from x,y to z,z, we can 
finally write down the above PDF of as 


dz 


t/2 


dcj) 


de 


Pt>,g{Xt = (Ttz,rfi) = \im I —f — / —pB{z,e,(j)) 


V—>-oo 


'-7r/2 


TT 


w —2N{l—cos(l>) —ipa-tzsin(l>—2pa-tQ(z){l—cosd)ipcrtzsm9 icpr —i9r 
A o e e e e e . 


( 20 ) 


where B{z,9,(l)) = f{z) + pat[2Ai{z)A2{z) cosfi + Al{z)e~^^ + Al{z)e'^^]. For large N, 
the major contribution of the integral over comes from the region around 0 = 0. 
Therefore, the 0 integral can be performed by expanding around 0 = 0 to make it 
a Gaussian integral (while extending the limits to ±cx)). Subsequently, one can also 
perform the Gaussian integral over This leads to the exact expression, 

Pt.g{Xt = atz.rfi) = - r (21) 

(Tt 271 

where 

B{z,9) = B{z,9,0) = f{z) pat[2Ai{z)A2{z) A‘l{z)e Al{z)e’'^y (22) 


In fact, using the integral representation of the modihed Bessel function of the hrst 
kind, for integer order n, 


d9 


In{x)= I [evidently, /_„(a;) = 4(x)] 


the above exact PDF given by (|21|) can be expressed in the closed form, 

'\/2Q(0y'+^ 


Pts^giXt = atZ.rfi) = — e 

(Tt 


-2patQ{z) 


X 


a/2Q(x) - X 

f{z) + 2patAi{z)A2{z) iJpcrt \/4Q^(z) - z^ 


+ pat A^(z) ^ ^ Ir+i (pat^/4:Q^{ 

^2Q{z)-z ^ 


z) — z^ 


+ pat Al{z) ^ Ir-I (pat^/4:Q^{ 

^^2Q[z) -\- z ^ 


z) — z^ 


(23) 


(24) 


3. Large deviation result for the two-tags displacement 


In this section, we obtain the large deviation form of (21) by evaluating the 9 integral 
using saddle point approximation. Note that in the expression of B{9,z) in (22), the 
second term is larger by 0 {pat) compared to the hrst term f{z), which comes from the 
process (i) where, in the non-interacting picture, the same particle happens to be the 


'th 


j"” and (j + r)^^ particles at the initial and final times respectively. So this process does 


not contribute at 0 {pat). 
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In this case we also scale the final with pat- The saddle point approximation of the 

1 


integral in (21) gives 


-Ptag(dft = atz, r = pad, t) 


(Zt ^j2Tig2{z) 


9i{z)e 


-patl{z) 


(26) 


where the different functions are explained below. 

The large deviation function (rate function) is given by 

I{z) = 2Q{z){l — cos 6*) — iz sin 6* + i6*l, 

where the saddle point 6 * is obtained using the condition 


d 


— \2Qiz){l — cos 9) — izsinO + i9l\ 
oO 


e=e* 


= 0 , 


(26) 


(27) 


which gives 




±1 +^/E + AQ\z)-z^ 

2Q{z)±z 


(28) 


Substituting 9* in (|26|) gives the large deviation function, explicitly in terms of as 

'i + ^/WTWW^ 


I{z) = 2Q{z) - + AQ‘^{z) - z^ + I In 


(29) 


2Q[z) + z 

The function gfiz) comes from evaluating the prefactor in ( [^ at the saddle point, 
gfiz) = {pat)~^ B{z,9*) 

= 2 Ai{z)A2 {z)A‘ l{z)e + 242(^)e*^ +0([pcri] (30) 


where are given in (28). 

The function g 2 {z) comes from performing the Gaussian integral around the saddle- 
point 9*, 

f)2 

g 2 {z) = ^—\2Q{z) {I - cos 9) - iz sin 9 + i9i\ = \/P + AQ‘^{z) - z'^. (31) 

o9^ e=e* 

The large deviation function has a minimum aX z = I and near this minimum, we 


get 




(32) 


2 2Q{1) 

Therefore, near the peak at Xt = r/p the PDF Ptag{Xt,r,t) has a Gaussian form, 
which describes the typical fluctuations. However, away from this central region, the 
Gaussian approximation breaks down, and one require the large deviation result (25) to 
describe for the atypical large fluctuations. In figure we compare both the Gaussian 
approximation and the large deviation result of the PDF with numerical simulation, and 
find that while the Gaussian approximation fits the data well near the central peak, the 
large deviation result agrees very well with the simulation data even beyond the central 
region. 
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Figure 6. (Color online) The (blue) points represent the simulation results for the 
PDF of the two-tag displacement in a one-dimensional system of hard-point particles 
for tag separations (a) r = 0, (b) r = 20, (c) r = 50, and (d) r = 100 respectively. The 
(magenta) thick dashed line corresponds to the analytic result in Eq. (25), while the 
(black) thin dashed line is for Gaussian distribution with the variance given by (37b I. 


4. Cumulants 


Now, we look at the cumulant generating function of the two-tag displacement X*. We 
dehne 

Z(A) = (33) 

such that the expansion of /i(A) in terms of the cumulants is given by 




pat ^ n\ 

’ n=l 


(34) 


Using the large deviation form of Ptag(Xt,r, f) given by (25), and then evaluating 


the integral over z; using the saddle point approximation, we have p{X) = A^* — I{z*) 
where z* is implicitly given by the equation A = dl{z*)/dz*. Using the expression of 
I{z) obtained above in terms of 6 * with the substitution = u we can express /i(A) 
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in the parametric form 


/i(A) = \z -\ -(z + /) + / In z/, 

/ 1 + z/^ ^ 

A = (l - z/"^) [1 + (z/ - 1) Ai{z )], 


V = 


-I + + AQ\z) - z^ 

2Q{z) - z 


(35a) 

(35b) 

(35c) 


The absence of p and a in (35) indicates that /i(A) does not depend on them. 


Therefore, from (34), it follows that 


W)c OC 


0-t 


p' 


,n—1 


(36) 


Note that the above equations (35), obtained through the saddle point calculation. 


gives the cumulants only at the most dominant order O(crt). To obtain the cumulants. 


we first expand the right hand side of (35c) about ^ = I and then invert the series to 


obtain 2 ; in terms of a series in z/ about z/ = 1. Therefore, the right hand side of (35b) 


can be expressed as a series in (z/ — 1). Next, by inverting (35b), we obtain z/ (and hence 


also z) in terms of a series in A about A = 0. Finally, from (35a), we express /i(A) as a 


series in A, and using the definition in (34), we obtain the first few cumulants (mean, 
variance, skewness, and kurtosis respectively) as 


a 


a 


p{Xl), = 2Q{1) 


a 


-1 ^2 


a. 


-1 


p^{Xi}, = 12F(l)Q(,l)-l, 

p’>{Xf), = 6 r8/(i)Q"(i) - Q(0 + 20F=(i)Q(i) - 21F{1 )], 


(37a) 

(37b) 

(37c) 

(37d) 


where I = r/{pat) and F{1) = f{w) dw. 

While the mean {Xt)c = r/p is exact to all order, the expressions for the other three 
cumulants are exact only at the leading order 0{at). The sub-dominant corrections can 


be computed using the exact expression of the PDF given by (21). For example, for the 
variance we get 


{XX = ^{2paiQ{l) + 




-1 


2f{l)Q{l) + 2F^{l)-- 


l 


Q\l)f"il) + fil) SF{l)Q{l) - - + /(/) 12 F\l) - 1 + 6 f{l)Q{l) 


0{[pat] 


(38) 


5. Velocity autocorrelations 

As a spin-off of our calculation, we show here that for the Hamiltonian model of 
elastically colliding particles, we can also compute the velocity auto-correlation function 
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Figure 7. The points are numerical simulation results for the velocity auto¬ 
correlation function computed from two tagged particles as a function of the tag 
separation at three different times, for the Hamiltonian model. The dashed line 
represent the analytical result to the leading order. 


(fo(0)t;r(^))- This can in fact be derived directly from the positional correlation function. 
We note that for the Hamiltonian case, we have Xi{t) = j^dt'vfit'). Hence it follows 


that 




1 d 

2 dt 


{xl(t)) + {xj(0)) - 2{x,(«)x„(0)) 


(39) 


The first two terms inside the square bracket are independent of time, hence they drop 
off on taking a time-derivative. For a Hamiltonian system we have {d/dt){xr{t)xo{0)) = 
{vr(t)xo{0)) = {vr{0)xQ{—t)). Taking another derivative, we get 


1 d^ 
2dt^ 


(llr)*) -lci(0)|“) = {tl,(()tlo(0)) . 


(40) 


Using Eq. (38), the fact that d'^Q/df 
we therefore get 


f(l), and at = vt for Hamiltonian dynamics. 


= (A-?) 


= —rf(l) + O 

pxt 


(41) 

(42) 


To leading order, this result can be obtained by a simple argument. Since the initial 
velocities are chosen independently for each particle, the contribution to the correlation 
function {nr(t)no(0)) is non-zero only when the velocity of the r—th particle at time t is 
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uo(0)ur(t)) ~ (5(r - pvt) v^) = ^ 

= - [ dfv-- 

pt J \ pt 


V -f [-) dv , 

V Vu/ 


hence hnally 




(43) 


(44) 


as in Eq. (42). In hgure we show a comparison of this analytic result with direct 
simulation results for the two-particle velocity autocorrelations in the equilibrium hard- 
particle gas. An exact expression for the velocity autocorrelation was obtained in [1] 
and involves a very lengthy calculation. This exact result can be recovered from using 
Eqs. ( 24|[4l ). However we see that the leading order expression is already quite accurate 
in describing the long time behavior. For the special case r = 0 and equilibrium initial 


conditions, using (38) in (41), we get 


(va(O)Vr(t)) ~ (pvt) “ (3^(0)/"(0) - /(O) + 6Q(0)P(0) 


(45) 


For the Gaussian distribution, f(x) = exp(—x^/2)/\'^, this gives (no(^)'?^o(0))/h^ ~ 
— {pvt)~^{2'K — 5)/(27r)^/^, a result hrst derived in [T], and also in (TB] using the present 
approach. 


6. Discussion 

We have considered a system of point particles moving on a one-dimensional line. The 
dynamics of individual particles is arbitrary and can be either stochastic or deterministic, 
and the only interaction between the particles is when they meet. The interaction 
dynamics is specihed by imposing that, we start with the non-interacting trajectories, 
and then interchange particle labels whenever trajectories cross — thus the ordering 
of particle labels is maintained at all times. This dynamics is quite natural for the 
deterministic so-called Jepsen gas and also for non-crossing Brownian walkers, and also 
seems natural for other stochastic processes. Using the fact that a mapping to non¬ 
interacting particles is available we have developed a formalism that seems to be very 
suited to computing reduced distribution functions and correlation functions in the 
interacting system. In particular, here we focus on computing the joint distribution of 
the positions of two tagged particles at different times. This is obtained exactly and 
from this we extract the large deviation function, various cumulants. For the case of 
Hamiltonian dynamics, we show that the two-pont velocity autocorrelation function can 
also be computed. We expect that the general strategy of our approach will be useful 
in the computation of more complicated correlations and distribution functions. 
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